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An optimal dynamical decoupling of a quantum system coupled to a noisy environment must 
take into account also the imperfections of the control pulses. We present a new formalism which 
describes, in a closed-form expression, the evolution of the system, including the spectral function 
of both the environment and control noise. We show that by measuring these spectral functions, 
our expression can be used to optimize the decoupling pulse sequence. We demonstrate this ap- 
proach with an ensemble of optically trapped uhra-cold Rubidium atoms, and use quantum process 
tomography to identify the effect of the environment and control noise. Our approach is applicable 
and important for any realistic implementation of quantum information processing. 
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Introduction - The paradigm of a two-level system (TLS) is central to quantum information (QI) , where it is applied 
as a quantum bit (qubit), the building block for information transfer, quantum memory or the computational gates. 
A Coupling to a noisy environment, which is inherent to any system, reduces the purity of the TLS, and thus limits 
its usefulness for any QI application. Several techniques have been developed to increase the quality of the quantum 
operation of the TLS, which is usually quantified by a measure called fidelity Pj. 

One of these techniques is dynamical decoupling (DD), where a pulsed control field is used to couple the two levels 
of the TLS, and thus reduce their coupling to the environment [5]. In QI, DD is used mainly to reduce the decay 
of the fidelity of a TLS, making it useful for longer times, as demonstrated in a large variety of systems [51-[TU]. It 
was theoretically and experimentally shown that the success of these schemes can be predicted using a measurable 
spectral function that describes the coupling of the system to the environment, sometimes referred to as the "bath 
spectrum" [inHTT] . In all these proposals, the higher the rate of the control pulses, the better is the decoupling 
from the environment. However, in any realistic implementation there are imperfections in the control field, hence 
dynamical decoupling sequences become increasingly ineffective as the number of pulses grow. 

In the quest for robust DD in the presence of pulse imperfections, it is usually assumed that errors in the control 
pulses are uncorrelated [TH]-[2D]. Such a "white" noise assumption is commonly used when estimating the fidelity of 
a gate, using the same gate many times but in a random manner, an approach known as "benchmarking" |21j . In 
contrast to this assumption, as we show below, the control pulses often have correlated errors leading to a non-flat 
spectral function, which must be taken into account. 

In this work, we study the combined effect of coupling to the environment and DD control field with colored noise 
on the TLS, and develop for it a closed-form expression, using the two corresponding spectral functions. We measure 
these spectral functions with an ensemble of ultra-cold optically trapped ^"^ Rh atoms, and then use them to predict 
the outcome of a generic DD scheme and its overall fidelity. Using quantum process tomography ^ Hj , we show and 
explain the effect of each of the spectral functions with any initial state. 

System subjected to realistic dynamical decoupling - We consider a general TLS model with energy fluctuations and 
a noisy control field, described by the effective Hamiltonian: 

H^^[uQ + 6{t)]cy, + ^[^l{t)e-''^°'o,^h.c] . (1) 

Here, wq is the transition frequency between the two states and f2(t)e~*"°* is a noisy (classical) external control field, 
which is used for the DD. The operators Gi are the Pauli matrices, written in the two level basis denoted by |2) and 
The noise in the control field enters through ri(<) = ilo(i) [1 + ?^c(^)], where ^o{t) is the desired, noiseless control. 
The frequency detuning noise, and the control noise, nc(t), are random functions of time, with zero average. Our 
analysis can be readily extended for control flelds having both multiplicative and additive noise, multi-axis pulses, 
and frequency and phase noise in the control [22] . 

The short time evolution of an initial state can be described by the reduced density matrix: p — \ {px<^x + Py<^y + Pz<^z + I), 
to second order in the noises |22]. In the interaction picture of Ho{t) = 1iax^o{t) /2 + 'haziOQ/2, and under the weak 
coupling assumption [Hj, the effective Hamiltonian, Hint{t) = ?i(5(i)^cos(/p dTl7o(r)/2) -|- hnc{t)^lo{t)^ , can be 
considered as a perturbation. The master equation for the density matrix operator is: 



dto 



H,nt{t),H„,t{t2)p{t) 



h.c. 



(2) 



where (...) stands for expectation value, after tracing out the environment. The two-term interaction Hamiltonian 
plugged into Eq. [2] gives rise to four terms, which are integrated over time T to flnd the short-time evolution of the 
density matrix. The outcome is the paper's theoretical main result: 
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(3) 

The three spectral overlap integrals in Eq. [S] determine the full evolution of the density matrix. They describe (in 
this order) the effect of the coupling to the environment, the noise in the control field and cross-correlation between 
the environment and the control field. The first term is similar to the spectral overlap integral of Refs. pTl [Ml [T7] . 
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The two bath spectral functions Gs{f) and Gc{f) (describe the correlation at different times of the environment 
and the noise of the control: 



Gsif) 
Gcif) 



— oo 
oo 
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-2-nifT 



{5{t) ■ S{t + r)) dT 



{nc{t) ■ nc{t + r)) dr , 



(4) 



The filter spectral functions Fs{f) and Fc{f) encapsulate the information regarding the modulation done by the 
control, during the time period T, and are written explicitly for sequences composed of or tt^^ pulses as: 



Fsif) 
Fcif) 



f dte-'^'"f^ ■ cos (^j^ no{T)dT^ 



(5) 



Similarly, Fsdf), Gsdf) and Gcs{f), describing the cross-correlations between the control noise, nc{t), and the 
environment noise, 6(t) are given in [22] , but are negligible in our experiment. 

Inverting the relation in Eq. [Sj in order to find the bath spectral functions from time evolution measurements of 
the density matrix is hard, when two or more overlap integrals are involved. However, it is possible to reduce the 
expressions to a single overlap integral by choosing wisely the initial state and DD sequence, essentially separating 
the problems of finding the two spectral functions. For example, to find the environment bath spectrum, Gs{f), we 
have used a random initial state and employed envelope spectroscopy to be insensitive to the control noise. With 
this choice, the evolution of the reduced density matrix, as given in Eq. 3, depends only on a single overlap integral 
[Hj. Using a filter function consisting of several discrete peaks, which samples the environment bath spectrum at 
these discrete frequencies, one can invert the spectral overlap integral solving a linear set of equations [HP or a single 
equation in the case of a single peak filter function |14) . 

Similarly, in order to measure the control noise spectral function, Gc(/), it is worthy to cancel the overlap integral 
of the environment. This is done by applying 7r-pulse sequence starting with an ensemble initialized to the state 
p(0) = Since it essentially keeps the system in states |1) and |2) that are insensitive to the pure dephasing 
environment noise, it reduces Eq. |3]to [22]: Ap(T) — ^ (py(T)) , with 



dfGdf)Fc{f) 



(6) 



By applying 7r/2-pulse followed by state detection, we can measure (^Py(T)'^, which is sensitive to the overlap integral, 
as explained. 

Measurement of the control noise spectrum - Our experimental setup is described elsewhere [23'. In short, ~ 2.5-10^ 
ultra-cold atoms are confined by an external optical potential created by two 1.06 /im crossed laser beams. 

The ensemble temperature is 1.7 fiK, and it has a central density of 2 • 10^^ cm~^. The two metastable states 
\1) = \F = l,m = —1) and |2) = \F — 2,m = +1) of the 5^5*1/2 manifold are chosen as the TLS. The energy difference 
between these states is, to the first order, magnetic insensitive, at the applied magnetic field of 3.2 G 24 . The control 
field is implemented using a two-photon MW-RF transition detuned by A ^ 110 kHz from the = 2, m = 0) level, 
taking into account all energy shifts (differential AC Stark shift, second order magnetic shifts, mean-field interaction 
and MW dressing) . We measure the state of the atoms using fluorescence detection scheme jl] . 

The main sources for the noises d{t) and ndt) are well understood in our system. The environment noise, 6{t), is 
due to the differential AC Stark shift of elastically colliding atoms in the optical dipole trap 4J. For each of the atoms, 
the environment is the atomic ensemble itself, randomizing the atomic trajectory following every elastic collision. The 
noise in the control is due to magnetic fluctuations. The magnetic noise enters through the single-photon detuning of 
the two-photon transition, A(t), which is first order magnetic field sensitive. This fluctuating detuning changes the 
effective Rabi frequency of the two-photon transition. 



n{t) = rjir22/2A(i) , 



(7) 



where fii and ft2, the single photon Rabi frequency of the MW and RF fields, are essentially noiseless in our system. 
Note, that since the two states of the TLS are magnetically insensitive, noise in the magnetic field affect only the 
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FIG. 1. The control noise spectrum, (a) The square (circle) data set is the measured (py) after 40 pulses of CPMG-4, with 
(without) injection of 4mG magnetic noise at 50 Hz. The solid lines are calculated using Eqs. |6]and[7) using the independently 
measured magnetic noise. Panels (b) and (c) : The filter function Fc{f) of a CPMG-4 sequence with pulse rate of 50 Hz and 
67 Hz, respectively. The arrows represent the 50 Hz component of the noise, which have significant overlap with the peaks of 
the filter function only for the sequence with a 67 Hz pulse rate (panel (c)). The inset shows an illustration of a CPMG-4 
control field at a 67 Hz pulse rate. The negative pulses are produced by a tt phase shift of the control field, (the duration of 
the pulses is not to scale). Points A, B, C and D are defined in the text. 



control field, hence the cross-correlation term in Eq. [3] is negligible. We expect to find a dominant contribution to 
the noise in 50 Hz and higher harmonics, arising from the electrical grid. 

In Fig. [l] (a) we plot the measured variance, (py), versus the pulse rate for 40 pulses of CPMG-4 DD sequence. 
The sequence CPMG-n, initially introduced by Carr, Purcell, Meiboom, and Gill is composed of equally 

spaced tt pulses vifith a phase alternating between tt and — tt after every n/2 pulses. We repeat the measurements with 
and without deliberately injecting a 50 Hz magnetic noise (by driving a current in a nearby coil, phase locked to the 
electrical grid) to farther increase the control noise. The measured spectrum is consistent with a single component 
at 50 Hz. Notice that since the filter function of CPMG-4 at pulse rate of 50 Hz has no peak at this frequency, 
there is no special feature around 50 Hz, as illustrated in Fig. [l] (b). This is in contrast to a pulse rate of 67 Hz 
which is overlapping a peak of the control noise, as shown in Fig. [l] (c). Using Eqs. [5][7]and a direct independent 
measurement of the 50 Hz magnetic noise, we depict the calculated control noise spectrum (in [l](a)), showing discrete 
peaks at frequencies of 200/(2m— 1) Hz, m = 1, 2, ... in excellent agreement with the measured noise spectrum without 
any fit parameter. The entire spectrum has a small bias which is due to imperfections in the state detection and 
uncorrelated (white) noise in the control pulses. The latter is measured separately using a higher number of pulses 
(100), to increase the sensitivity. Finally, a DC control noise of 3 • 10^^, is measured using a CPMG DD sequence 
(with no phase alternation), whose filter function has a prominent component at DC. 

By choosing one of the peaks in the spectrum of[l](a) and repeatedly measuring its noise component, we were able 
to reduce it by about 50 percent by injecting 50 Hz component to the nearby coil and searching for the phase and 
amplitude which minimize the peak. 

DD Sequence engineering - The usefulness of Eq. [3] stems from its ability to predict the performance of any DD 
sequences, given that the two spectral functions, Gs{f) and Gdf), are known. For the purpose of optimizing the 
DD we quantify its success with the fidelity, defined by J" = Tr{p{t)p{0)}, as it includes both the effects of pulse 
imperfections and the coupling to the environment. The short time fidelity (neglecting the last term in Eq. [3|, can 
be written as [22] : 

2 ^ 2 ^oo 

^ J — oo 
2 I 2 .oo 

/ dfGM)Fc{f) • (8) 

^ J — oo 

It is helpful to plot the decay rate of the fidelity as a function of the pulse sequence parameters. In such a plot, it is 
easy to graphically identify the region in parameter space where the performance of the DD sequence is optimized. In 
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FIG. 2. A two-dimensional map of the predicted logarithm of fidelity decay rate (in units of sec~^) presented in the variables 
pulse rate, fno and phase cycle n of the CPMG-n pulse sequence. The shade (color online) represent the decay rate of the 
fidelity in logarithmic scale (black represents no decay). The marked zones of best expected performance of the DD sequence 
are identified by eye inspection. Points A, B, C and D are defined in the text. The dashed line correspond to n = 4 (i.e. 
CPMG-4 sequence). 



the case of CPMG-n, the natural choice of parametrization is the pulse rate, fjjo and the parameter n. An example 
for such a calculated map, based on the two measured spectral bath functions in our system is presented in Fig. [2] 
for the worst case fidelity (taking an initial state py — 1). The regions with the highest fidelity are clearly visible. 

Although the map exhibits some complex features, the central ones can be qualitatively understood. At low control 
pulse rate, the fidelity decay rate follows a Lorentzian, reflecting the Poisson statistics of the cold atomic collisions 
[3]. For higher pulse frequencies, there is a reduction in fidelity owing to white noise in the control field arising from 
the large number of pulses. Large-n cycles are also less successful, since their filter function has a large spectral 
component in DC, sampling the slow drifts in strength of the control field. The control noise, originating mainly from 
the 50 Hz magnetic noise, produces a dominant feature appearing as strips on the map. The points, A-D correspond 
to the points marked in Fig. [T] 

Coherence of an arbitrary initial state - Process tomography [J is a technique used to characterize the TLS state 
after being manipulated by the control (referred here by process), for any initial state. A great advantage of the 
formalism presented in Eq. |3]is that it predicts the entire 3-dimensional effect of the process on the system, which is 
simply visualized as deformed sphere in the Bloch representation. 

For process tomography we repeat the process with four initial states, p(0) = ^ (S -I- 1), where E is the Pauli matrix 
Cz, —f^zi or ay. For each initial state, the final state is measured by applying 6 different control pulses, followed 
by a state detection. For a linear process, this information is sufficient to construct the process matrix [Ull]. 

The results of the process tomography measurements are shown in Fig. [3j for two DD processes: CPMG-4 at 50 Hz 
and 67Hz pulse rates, corresponding to points B and C in Fig. [2j The two processes, although close in frequency, differ 
significantly, in agreement with our measured bath spectral functions. For the 50 Hz process, there is no dominant 
control noise (as explained before), hence the decay is mostly due to the coupling to the environment. The decay rate 
of the z axis is the slowest, limited by a Tl process (not included in our model, measured to be ^ 5sec). The decay 
of the other two axes is similar, which is expected from Eq. [3j since the coefficients px and py appear symmetrically 
in the terms describing the coupling to the environment. In contrast, in the 67 Hz process, the decay in the y and z 
axes is faster since it is also affected by the noise in the control (involving the coefficients py and py). 

Conclusions - The formalism developed here together with the one described in '14] gives a recipe for designing a 
DD sequence: First measure the spectral function defining the coupling to environment. Then measure the spectral 
function that characterizes the noise in the control field. Choose a general DD sequence parameterized by few 
parameters. Use the overlap integrals to calculate the performance map as a function of these parameters. Choose 
high fidelity regions for the DD sequences. This framework can be extended also to sequences combining and ny 
pulses. The filter functions, however, don't have an analytical expression, in this case. Although the the optimal 
pulse is system dependent. The most efficient way to find it is general, as we have shown here. 

We acknowledge the financial support of MINERVA, ISF and DIP. 
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FIG. 3. Process tomography of CPMG-4 DD scheme in the presence of a noisy control field. A Color code is used to help 
determine the outcome with any initial state. Top: Measured Bloch sphere after 50 pulses at a rate of 50 Hz with 4 mG 
magnetic noise injection. Top and side view on the right and left respectively. Bottom: The same with 67 Hz DD. The different 
axes suffer from different decay rates in each case as expected from Eq. |3] 



[1] M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information (Cambridge University Press, 2000). 
[2] M. H. Levitt, Spin Dynamics: Basics of Nuclear Magnetic Resonance (Wiley-Blackwell, 2008). 
[3] S. Kotler, et al, Nature 473, 61 (2011). 

[4] Y. Sagi, I. Almog, and N. Davidson, Phys. Rev. Lett. 105, 053201 (2010). 
[5] L. Viola, S. Lloyd, and E. Knill, Phys. Rev. Lett. 83, 4888 (1999). 
[6] L. Viola, E. Knill, and S. Lloyd, Phys. Rev. Lett. 82, 2417 (1999). 
[7] C. Search and P. R. Herman, Phys. Rev. Lett. 85, 2272 (2000). 
[8] K. Khodjasteh and D. A. Lidar, Phys. Rev. Lett. 95, 180501 (2005). 
[9] G. S. Uhrig, Physical Review Letters 98, 100504 (2007). 
[10] J. Bylander, et al, Nat Phys 7, 565 (2011). 

[11] A. G. Kofman and G. Kurizki, Phys. Rev. Lett. 87, 270405 (2001). 

[12] G. Gordon, N. Erez, and G. Kurizki, Journal of Physics B: Atomic, Molecular and Optical Physics 40, S75 (2007). 

[13] H. Uys, M. J. Biercuk, and J. J. Bollinger, Phys. Rev. Lett. 103, 040501 (2009). 

[14] L Almog, et al, Journal of Physics B: Atomic, Molecular and Optical Physics 44, 154006 (2011). 

[15] G. A. Alvarez and D. Suter, Phys. Rev. Lett. 107, 230501 (2011). 

[16] T. Green, H. Uys, and M. J. Biercuk, Phys. Rev. Lett. 109, 020501 (2012). 

[17] J. Clausen, G. Bensky, and G. Kurizki, Phys. Rev. Lett. 104, 040401 (2010). 

[18] L. Viola and E. Knill, Phys. Rev. Lett. 90, 037901 (2003). 

[19] A. M. Tyryshkin, et al, (2010), arXiv:1011.1903 [quant-ph]. 

[20] A. M. Souza, G. A. Alvarez, and D. Suter, Phys. Rev. Lett. 106, 240501 (2011). 

[21] E. Magesan, et al., Phys. Rev. Lett. 109, 080505 (2012). 

[22] For more details see the supplementary material section. 

[23] Y. Sagi, L Almog, and N. Davidson, Phys. Rev. Lett. 105, 093001 (2010). 

[24] D. M. Harber, et al, Phys. Rev. A 66, 053616 (2002). 

[25] H. Y. Carr and E. M. Purcell, Phys. Rev. 94, 630 (1954). 

[26] S. Meiboom and D. GiU, Rev. Sci. Instrum. 29, 688 (1958). 



Supplementary material for 
Optimal dynamical decoupling in the presence of colored control noise 



Ido Almog^, Gil Loewenthal^, Jonathan Coslovsky-'^ , Yoav Sagi^, and Nir Davidson'^ 
^Department of Physics of Complex Systems, Weizmann Institute of Science, Rehovot 76100, Israel and 
^ JILA, National Institute of Standards and Technology and University of Colorado, 
Department of Physics, University of Colorado, Boulder, Colorado 80309-0440, USA 



2 



EQUATION OF MOTION FOR SYSTEM DENSITY OPERATOR 
Mathematical derivation 

We are interested in the short time evolution of the system density operator determined by the Hamiltonian: 

H=^[uJo + S{t)] a, + ^ [0(t)e-^"°V, + h.c] . 



(1) 



It is convenient to change to the interaction picture of Ho{t) = hax^o{t)/2 + Tiaz^o/'^- In this picture, when the 
change in the system density operator during the bath correlation time is negligible, the effective Hamiltonian, Hint, 
can be considered as a perturbation. This assumption, called the weak coupling assumption, is naturally valid in QI 
since the quantum gate or memory is not useful when the decay of coherence is not small comparing to unity. In this 
case, the system's density obeys the following master-equation (to second order in the noises): 



m 



(in) 



h)'Jo 



dto 



Hint{t), Hint{t2), p{t) 



(2) 



where (...) stands for expectation value, coming out from tracing out the environment. 

The interaction hamiltonian can be written explicitly by Hint{t) = h5{t)Ss{t) + Jinc{t)Sc{t), with an operator 
responsible for the interaction with the environment: 



^cosm) + y .^m(a?(i)) 



and an operator responsible for the noise in the control 

Sc{t) 



i(Tx X —i<7x 



(3) 



(4) 



where the angle 'd{t) = Jg dTQo{T)/2 is sometimes called the action. Plugging this into Eq. 2 we get the equation for 
the density operator: 



-f 

Jo 

f 

Jo 



dt2{6{t)S{t2)) Ssit), Ss{t2),m 



dt2 {nc{t)5{t2)) Sc{t), Ss{t2), pit) 



f 

Jo 



'Sc{t2), 


m] 


4(i2), 


m] 



(5) 



Three mechanisms responsible for the change of the density matrix can be identified from the four terms in the 
expression: The first term describes the effect of the environment, the second describes the effect of the noise in the 
control field and the last two are mixed effect resulting from correlation between the control noise and the environment 
coupling. Fortunately, each term in Eq. 5 can be expressed as an overlap between two function in the spectral domain. 
For example the last term can be written: 

1 



- / dt2 {5{t)n,{t2)) Ssit), S,it2)m 
Jo 



dt2 {Sit)n,it2)) -aMt2) [PyiO) sin (-^(f)) + p,(0) cos (-^(f))] .(6) 



We integrate to get the expression for density matrix change contributed by this term: 



ApsciT)= [ dti ( dt2^Sciti-t2)-laMt2)pziO)cosi^iti)) 
Jo Jo ^ 



(7) 



Note that to get this we ignored the sin(...) term since we constrain ourselves to sequences consisting on tTj, pulses 
only. Here we used the definition for the correlator which is a function of the time difference only: $5c(ii — ^2) = 

((5(ti)nc(f2)) 0(^1 ~ t-i), with heaviside fimction 9(fi — t'l)-, looking forward to the next step where we replace the 
upper integral limit. Identifying a convolution integral between the functions $5c(i) and fio(i)i with the notation 
S(f) = cos i-dit)) e(i)e(T - t), we write: 



Ap5c(T) = \oxpM ^ c!iiS(ti) ($5c(ti) * ^oitx)) 



(8) 
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Using the Fourier transform of the correlator (defined with a 2 factor to compensate for that the transformation of 
only the positive times) Gsdf) = 2 dTe~'^^^^'^<^sc{T) and the symbols ~ for Fourier of other functions, the integral 
(without the coefficient) becomes 



Ap5c(T) = -j^ dhE{h)^s,ih) * noih} = - y dhE{h) J dhe^^'f^'^GUfi) ' ^^o(/i) 

1 /"CSO /"OO 1 fOO 

= dt d/2e2-/^*S(/2) / d/ie2-/^*G,e(/i)Oo(/i) = - / dfE*{f)no{f)Gsc{f) • (9) 

J — oo J —oo J — oo J — oo 

With the notation Fsc{f) = S*(/)17o(/), the density matrix change resulting from this contribution: 

Apsc{T) = Ict^pM r dfGsc{f)Fsc{f) ■ (10) 

^ J — oo 

Repeating this procedure taking the three other terms in Eq. 5, obtain the other overlap terms. The first two 
terms, describing the uncorrclatcd effect of the coupling to the environment and the noise in the control, are written 
explicitly in the Letter. The forth term is similar to Eq. 10 and can be written as: 

1 f°° 

Apcs{T) = -a^pM / dfGcs{f)Fso{f) , (11) 

^ J — oo 

with Gcsif) defined the same as Gsdf) to exchanging the function 6{t) and ndt). 

Generalization for other cases 

Frequency (and phase) noise of the control arc treated the same as the energy fluctuation S{t). To see this, one 
has to change to the interaction picture which rotates with the frequency of the control local oscillator, including its 
classical frequency noise. In this interaction picture, the noise in the control is an additional ct^ noise (up to a minus 

sign) added to S{t). 

Control pulses involving the two axes x and y should be considered carefully. Although changing the axis is 
equivalent to a phase shift, it cannot be regarded as energy fluctuations like is done with the frequency noise of the 
control, since the phase shift cannot be considered small, breaking the weak coupling assumption. Instead, one should 
carry on the similar derivation derivation as was done above. Applying this procedure, the commutators in Eq. 5 
become somewhat more complicated to calculate since there are no close expression for Ss{t) and Sc{t), since the pulse 
action now includes the two non-commuting operators, and (Ty. Nevertheless, Ss{t) and Sc{t) can be calculated 
directly, by evaluating the effect the pulses, sequentially. Since between each pulse the evolution can be calculated 
analytically, the resulting fllter function can be written as a as long table of function as the number of pulses, making 
it still possible to calculate the overlap integrals. In the case of continues control, the sine term we disregarded in the 
treatment should be taken into account like is done in [1]. 

The model for the noise in the control, modeled by fl{t) = flo{t) [1 + nc{t)], can be considered as general l^*-order 
treatment of the noise. In this model the noise is gated (multiplied) by the function flo{t)- However, one can consider 
also the case of additive noise which is not gated by this function: 

n{t)=no{t)[l + ne{t)] + b{t) (12) 

This noisy effect is the case if the control is derived from an oscillator which is controlled by an analogical gate which 
is noisy. However the effect of this noise, b{t) has the same efi^ect as nc{t), but gated with a function which is constant 
during the entire experiment time, T. To regard this in our treatment the filter function must be replaced by: 



Pcif) 



L 



T 

dte-'^'''f^ 



(13) 



which is just the original filter function plugging = 1. In the case where the two kind of noises appear simultane- 
ously, one should add an additional overlap integral, (or two additional overlap integrals in the case the noises b{t) 
and nc{t) are correlated.) 
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SHOT-TO-SHOT FLUCTUATIONS AND THE SYSTEM DENSITY OPERATOR 

The evolution of the density matrix, described in our framework, is a sufficient description also for the fluctuation 
of the experimental results. For example, a strong noise in the control field can completely randomize the final state 
of the system. Consequently, the parameters px, Py, and pz, which parameterize the density operator, have to be zero, 
for this randomized state. Similarly, measuring the shot-to-shot fluctuation, one can calculate the average parameters: 
px, Py, and pz- It is important to note that unlike the noise, the random effect of the control noise docs not average 
out on the ensemble since all atoms are subjected to the same noise in a single realization. Therefore, more than a 
single experimental measurement has to be done to know the parameters p^, py, and p^- 

The density operator fluctuations, arc related to the average change in the density operator. This relation is 
captured by a geometrical interpretation of summing random vectors: Taking the sum of random vectors pointing in 
random direction, one gets a vector which length is smaller than the sum of all vector length. Describing the density 
operator by a vector and applying this argument, one can derive the relation between its mean value and fluctuations. 
However, we here derive it similarly to the previous section, starting from the evolution of the density operator. The 
master equation for the density operator up to the 1^* order in the control noise is: 

p'e(t) = ^ [Hintit), m] ■ (14) 

Plugging the interaction Hamiltonian, integrating over the short time T and taking Trace to find the y-component of 
the vector, we obtain: 

ApyiT) = Tr{Ap,{T)ay} = Trj^ ^ dtn,(t) [3^(1), p(t)] tTj^j 

= / dtnc{t)pMt) . (15) 
Jo 

Taking the expectation value of Apy(T) yields zero as the noise nc(t) averages to zero. However the fluctuations in 
this component is readily calculated by: 



(Apy(r))^^ = dhnc{h)pzflo{h)- dt2nc{t2)Pz^o{t2)j , (16) 
which in the frequency domain is written: 

(Ap,(T))'^ = pI dfG^imif) , (17) 



where the density operator change, Apy{T), can be replaced by the component, Py{T) itself in this expression, since 
the initial density operator is not fluctuating. 

SHORT-TIMES FIDELITY 

The fidelity T = Tr{p{t)p{0)} is derived for short times for an initial pure state. For pure dephasing, it is written 
by: 

T = Tr{p{t)p{0)}=Tr{{p{0) + Ap{T))-p{0)} . (18) 
Since p(0) is a pure state, the trace of p^(0) is 1 and therefore: 
jr=l+Tr{Ap(T).p(0)} 

P.<rx + PyCJy T 1^ ^^G,(/)F,(/)) \ {pxCJx + PyOy + p.CT, + 1) | 



= l + Trl [-' 



2T 



= 1 - -Tr{al}pl— j ^ dfGsimif) - -Tr{al}pl— j ^ dfGs{f)Fs{f) 



5 



= 1 - ^ / dfGsimif) - -f / dfGsif)Fsif) . (19) 

Doing the same arithmetics with the noise in the control, one gets the fidehty of the system subjected to the two 
noises. It is written as: 

J- = 1 - ^^^Ti?5(T) - '^^^^TR^iT) (20) 
with the definitions for the decay rate functions: 

RsiT) = ^ dfGsif)Fsif) , (21) 

and 

i?c(T) = ^|_"rf/Ge(/)Fc(/) . (22) 

The definitions of RsiT) and RdT) are useful for since the functions are almost time independent for periodic pulses. 
This is due to the scaling of the filter function canceling the 1/T dependence. Hence the fidelity decays in a constant 
rate. It also helpful to note that the integral of filter function Fgif) is just T and theretofore, doesn't depend on the 
pulse frequency, but only on the total experiment time. However, the filter function F^f) scales linearly with the 
pulse rate. Hence, the decay rate of the fidelity is large for spectral features at higher frequencies which is consistent 
with the experimental results. 
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